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Abstract 

We present updated results on the search for a neutrino signal from the core of the Earth and of the Sun induced 
by Weakly Interacting Massive Particles (WIMPs). In this paper we concentrate on neutralinos as WIMP 
candidates. The 971 and 642 events used respectively for the search from the Sun and from the Earth are 
compatible with the background of atmospheric neutrinos. Consequently we calculate flux limits for various 
search cones around these sources. Limits as a function of the neutralino mass are given and compared to the 
supersymmetric (SUSY) models. 

1 WIMPs from the Earth and the Sun: 

We describe an indirect method to look for non-baryonic dark matter WIMPs. The best WIMP candidate is 
the lightest SUSY particle (LSP) which in the Minimal Supersymmetric Standard Model (MSSM) is expected 
to be stable if R-parity is conserved and hence should be present in the Universe as a cosmological relic from 
the Big Bang. The linear combination of Higgsinos and gauginos, the neutralino x, is currently considered 
the best candidate for cold dark matter since its couplings and mass range naturally give the relic density 
required to explain halo dark matter. Neutralinos are described by 3 parameters in the MSSM assuming the 
GUT relation between gaugino masses: one of the gaugino masses, the Higgsino mass parameter \i and the 
ratio of the Higgs doublet vacuum expectation values tan f3. Moreover, if universality is assumed, the MSSM 
phenomenology is described by these parameters, the universal trilinear scalar coupling A, the degenerate 
scalar mass ttiq and the mass of the pseudoscalar neutral Higgs rriA- Experimental searches at LEP set a 
lower limit on m x at ~ 32 GeV and suggest an upper limit at ~ 600 GeV if one requires that the neutralino 
Q\ | cosmological abundance £l x h? < 0.3 (Ellis, 1999). This upper limit is not yet achievable by the forthcoming 
LEP and Tevatron runs. Direct and indirect searches at underground detectors explore SUSY parameters and 
are complementary to the future LHC measurements. Direct searches look for a signature of a direct scattering 
of a WIMP from a nucleus in the detector. The DAMA experiment (~ 100 Kg Nal(Tl)) sees an indication 
of an annual modulation of the rate which could be due to the Earth's motion around the Sun and the change 
of the Earth's velocity relative to the incident WIMP. The 19511 kg day data favor at 99.6% c.l. the presence 
of an annual modulation signal which, if interpreted in terms of WIMPs, implies a mass of my/ = (591^) 
GeV (Bernabei et al., 1999). This indication should be checked using different techniques, such as the indirect 
detection of trapped WIMPs inside the core of the Earth and of the Sun. The signature would be an excess of 
neutrino events resulting from WIMP-WIMP annihilations around the direction of the vertical of the apparatus 
and of the Sun beyond the known atmospheric v background (Jungman, Kamionkowski & Griest, 1996). 
MACRO measures neutrinos indirectly as upward-going muons and has presented results of the WIMP search 
in Ambrosio et al., 1998a, to which we refer for details. We update this search including the data collected 
during Mar. 98-Feb. 99. 

2 MACRO Updated Results on WIMPs: 

The MACRO detector at the Gran Sasso Laboratories, with overall dimensions of 12 x 76.6 x 9 m 3 , detects 
upward-going muons through the time-of-flight measurement using 600 tons of liquid scintillator inside 12 
m long boxes (time resolution ~ 500 psec). A system of around 20, 000 m 2 of streamer tubes reconstructs 
tracks with angular resolution < 1°. The lower part of the apparatus is filled with rock absorber setting a 1 
GeV threshold for vertical /xs. The upward-going muon measurement relative to the construction period of 
MACRO (Mar. 89 - Apr. 94: 1.38 yr of running of 1/6 of the lower apparatus and 0.41 yr of the lower detector, 



inefficiencies included) is described in Ahlen et al., 1995. Since then, MACRO is in its full configuration (3.93 
yr, inefficiencies included) and further results are in Ambrosio et al., 1998a and Ambrosio et al., 1998b. 

For the WIMP search for the Earth we use the sample of 642 throughgoing upward muons selected with 
the requirement that the track crosses at least 200 g/cm 2 in the MACRO rock absorber, which reduces the 
background due to soft 7rs produced at large angles by downward-going ^s to ~ 1% (Ambrosio et al, 1998c). 
Releasing this cut, we use 971 upward-going ^s for the search for the Sun because background rejection is not 
so critical for moving sources and the increase in exposure offsets the slight increase in background. 



For the Earth, the expected background due to interactions of atmospheric z^s in the rock below 
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anu ine muon en- Table 1; observed and atmosp heric ^-induced background and 90% c.l. /j flux limits as a function of half- 
ergy loss as in Loh- cone angles around the Earth core and the Sun. For the Earth, the expected background events are multiplied 
mann et al., 1985. by the ratio of observed to expected events outside each cone. The Earth results are for the no oscillation 
We estimate a to- scenar i°- The average exposure for the Earth is 3272 m 2 yr and for the Sun 1 1 16 m 2 yr. 



tal uncertainty in the calculation of upward-going muon fluxes of 17%. For the Earth we have even considered 
a — > v T oscillation scenario with parameters Am 2 = 0.0025 eV 2 , sin 2 20 = 1 as suggested by the flux 
measurement reported at this conference (Ronga et al., 1999). For the Earth search, the expected number of 
atmospheric induced events in the no oscillation scenario is (including the contribution due to ^-interactions 
inside the bottom part of the apparatus which are selected as throughgoing muons) 835 ± 142 and in the os- 
cillation scenario 581 ± 99. 

For the Sun we have compared the 971 measured events with a different simulation with respect to that used 
for the Earth. This is obtained by mixing randomly the local coordinates of measured upward-going events 
and times gathered during the entire data-taking. This method takes into account the contribution of events 
produced by internal ^-interactions in the MACRO absorber. In Fig. 1(a) we show the angular distributions 
of the measured and expected events from atmospheric neutrinos for the Earth and the Sun. For the Sun this 
distribution depends strongly on time due to the motion of the Sun. In Fig. 1(b) we show the muon flux limits 
for 10 search cones from 3° to 10°. In the case of the Earth, the expectation from atmospheric us in the region 
of interest for the signal is larger than the data; we then evaluate flux limits multiplying the expected number 
of events by the ratio of the data to the expectation outside the cone where we look for the signal. This nor- 
malization is motivated by the high uncertainty in the normalization of upward-going muon flux calculations, 
whereas the shape error in the flux distribution is a few percent only. Moreover, since the number of detected 
events is less than the normalized expected events, we set conservative flux limits assuming that the number 
of measured events equals the number of expected ones (Caso et al., 1998). With this normalization, Earth 
limits considering ^-oscillations agree with the ones in the case of no oscillations within 7%. On the other 
hand, for the Sun, having used data to evaluate the expected numbers from atmospheric z^s, oscillation effects 
are automatically included in the given limits. In Table 1 we show measured and expected events and flux 
limits for some of the cones calculated assuming a minimum /i energy of 1.5 GeV and 2 GeV for the Earth 
and the Sun, respectively. The minimum energy for the Sun is higher because tracks pointing toward it are 



more slanted than vertical tracks pointing to the core of the Earth and hence cross a larger amount of MACRO 
absorber. The average exposures (live-time times detector area in the direction of the expected signal from 
the source of WIMP annihilation) for cones between 3° and 30° is 3272 m 2 yr for the Earth and 1116 m 2 yr 
for the Sun. We estimate a maximum error of 5% on flux limits assuming these minimum energies for flux 
limit calculation with respect to a calculation which takes into account the dependence of the acceptance of 
the apparatus and of the neutrino fluxes from x ~ X annihilation. 
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Figure 1: (a) On the left: zenith angular distribution of measured and expected events for the Earth search. 
Solid and Dotted lines: expectations for atmospheric us for no-oscillations and — > v T oscillations with 
Am 2 = 0.0025 eV 2 ; black circles: data. The expected distributions are multiplied by the ratio R of the mea- 
sured events over the expected ones outside the largest window 30° (R = 0.83 for the solid line and 1.15 for 
the dashed line). On the right: distribution of the cosine of the angle from the Sun direction, (b) /j, flux limits 
(90% c.l.) vs. the half-width of the cone from the source direction. 



We calculate flux limits assuming the neutralino is a WIMP candidate. We estimate search cones con- 
sidering the angular distribution of upward-going muons around the direction of neutrinos from neutralino 
annihilation which mainly depends on the neutralino mass using Bottino et al., 1995 flux calculations. We 
choose to calculate flux limits in those cones which include 90% of the expected signal. We estimate that 
extreme assumptions on the annihilation channels (the final states are fermion pairs, gauge/Higgs bosons) 
produce a maximum variation of flux limits of about 17%. The Bottino models are shown as dots (model 
parameters are described in detail in Ambrosio et al., 1998b) and circles in Fig. 2 as a function of neutralino 
mass for the Earth and the Sun and a minimum muon energy of 1 GeV. The solid lines represent MACRO 
flux limits corrected for the 90% collection efficiency of the cones and including correction factors to translate 
limits in Table 1 to the 1 GeV threshold used in the calculation. These factors are maximum for low neutralino 
masses (1% for the Earth and 10% for the Sun at m x = 60 GeV). As shown in detail in Ambrosio et al., 
1998b, MACRO'S experimental limit from the Earth rules out a considerable number of SUSY configurations 
based on the interpretation of the DAMA/Nal data (Bottino et al., 1999), even assuming neutrino oscillations 
of the atmospheric v background. It should anyway be considered that even if the MACRO flux limits vary 



within 7% if the oscillation or no oscillation hypothesis is assumed, — v T oscillations with the already 
quoted parameters could lower neutrino flux calculations from neutralino annihilation by at most a factor of 
two (Fornengo, 1999). In the case of the Sun, MACRO has less overlap in sensitivity for low neutralino masses 
with direct searches because the indirect search is more sensitive to spin-dependent scattering. 
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Figure 2: On the left: dots are upward-going \i fluxes vs m x for E™ = 1 GeV from the Earth for various 
model parameters (Bottino et al., 1995); solid line: MACRO flux limit (90% c.L). Open circles: models 
excluded by direct measurements (Bernabei et al., 1996). On the right: the same as on the left for the Sun. 
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